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We investigate the role of the initial condition used for the hydrodynamic evolution of the system 
formed in ultra-relativistic heavy-ion collisions and find that an appropriate choice motivated by the 
models of early-stage dynamics, specifically a simple two-dimensional Gaussian profile, leads to a 
uniform description of soft observables measured in the Relativistic Heavy-Ion Collider (RHIC). In 
particular, the transverse-momentum spectra, the elliptic-flow, and the Hanbury-Brown-Twiss cor- 
relation radii, including the ratio Rout/Rside as well as the dependence of the radii on the azimuthal 
angle (azHBT), are properly described. We use the perfect-fluid hydrodynamics with a realistic 
equation of state based on lattice calculations and the hadronic gas at high and low temperatures, 
respectively. We also show that the inclusion of the partonic free-streaming in the early stage allows 
to delay the start of the hydrodynamical description to comfortable times of the order of 1 fm/c. 
Free streaming broadens the initial energy-density profile, but generates the initial transverse and 
elliptic fiow. The data may be described equally well when the hydrodynamics is started early, or 
with a delay due to partonic free-streaming. 
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The notorious difhculties in simultaneous description 
of various features of the soft hadron production in the 
nucleus-nucleus collisions at RHIC, which occur within 
the standard approach consisting of partonic, hydrody- 
namic, and hadronic stages, are well known In partic- 
ular, one of the so called RHIC puzzles refers 
to problems in reconciling the large value of the ellip- 
tic flow coefficient, V2, with the Hanbury-Brown-Twiss 
(HBT) interferometry in numerous approaches including 
hydrodynamics [1, [a, B 111- In niost existing analy- 
ses, the large value of V2 prefers long evolution times of 
the system, while the HBT radii indicate that this time 
should be short. Typically, the hydrodynamic evolution 
is initiated from an initial profile generated by Glauber- 
like models, with the initial temperature serving as a free 
parameter. 

In this Letter we show that the choice of the initial 
condition for hydrodynamics is a very important element 
for the proper description of the data. Ideally, the ini- 
tial density and flow profiles should be provided by the 
early partonic dynamics, for instance the Color Glass 
Condensate (CGC) [l^, [HI. In practice, however, the 
theory of the partonic stage carries some uncertainty in 
its parameters, moreover, there may exist other effects, 
see e.g. [H, which influence the early dynamics. It 
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is then practical to use some simple parameterization of 
the initial profile. Here we investigate boost-invariant 
systems, which approximate well the RHIC collisions at 
mid-rapidity. We take the following Gaussian parame- 
terization for the initial density profile of the system in 
the transverse plane {xo,yo) at the initial proper time 
To = 0.25 fm, 

n(.o,yo) = exp(-||-f). (1) 

The values of the a and b width parameters depend on 
the centrality and are obtained by matching to the results 
for (x^) and (y^) from GLISSANDO f]3 . w hich implements 
the shape fluctuations of the system [l5| . The values for 
centrality classes used in this work are collected in Ta- 
ble m The profile ([T]) determines the energy-density pro- 
file, from which the temperature profile is obtained [la |. 
The initial central temperature, which may depend on 
the centrality, is denoted by Tj and is a free parameter 
of our approach. 

The hydrodynamics used is inviscid, baryon-free, and 
boost-invariant. The equation of state is taken to be as 
realistic, as possible. We use the lattice QCD simula- 
tions of Ref. [13] at high temperatures, T > 170 MeV, 



TABLE I: Shape parameters for various centrality classes. 



c[%] 


0-5 


0-20 


20-30 


20-40 


a [fm] 


2.65 


2.41 


1.94 


1.78 


b [fm] 


2.90 


2.78 


2.52 


2.45 



2 



1000. 




0.2 0.25 0.3 0.35 0.4 0.45 0.5 



c =20-30 % 

Ti =460MeV 

T, =145 MeV 

To = 0.25 fm 

T = 0.25 fm 



= 20-30 % 
= 305 MeV 
= 145 MeV 
= 0.25 fm 
= 1 .0 f m 





K + K 




TT + K 




P 




P 


c 


= 20-40 % 


c 


= 20-40 % 


Ti 


= 460 MeV 




= 305 MeV 


T, 


= 145 MeV 


T, 


= 145 MeV 


To 


= 0.25 fm 


To 


= 0.25fm . 


T 


= 0.25 fm 


T 


= 1.0 fm>0 




0.01 

0.001 
0.25 

0.2 
0.15 

0.1 
0.05 



Pt [GeV] 

FIG. 1: (Color online) The transverse-momentum spectra 
of pions, kaons and protons for the centrality class c=0-5% 
(the upper panel) and c=20-30% (the middle panel), and the 
transverse-momentum dependence of the elliptic flow coeffi- 
cient 112 for the centrality class c=20-40% (the lower panel). 
The darker (lighter) lines/bands describe the model results for 
the case without (with) free-streaming. The PIfENIX data 
are taken from Refs. [itI. [isl]. 



the hadronic gas at T < 170 MeV, and a smooth interpo- 
lation in the vicinity of 170 MeV, as described in Ref. 
According to recent knowledge, no first-order transition 
is implemented, but a smooth cross-over. The hydrody- 
namic equations are solved by the method of Ref. ^21j. 
The accuracy of the method is tested with the entropy 
conservation, satisfied at the relative level of 10~^ or bet- 
ter. At the temperature Tf — 145 MeV (model param- 
eter) the system freezes and hadrons (stable and reso- 
nances) are generated according to the Cooper-Frye for- 
malism. Possible elastic rescattering processes among 
these hadrons are neglected, thus they stream freely, with 
the resonances decaying on the way. This stage is simu- 
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FIG. 2: (Color online) The transverse-momentum dependence 
of the HBT radii -Rsidc (a), -Rout (b), iiiong (c), and the ratio 
-Rout / -Rsido (d) for central collisions. The darker (lighter) lines 
describe the results without (with) free-streaming. The STAR 
data are from Ref. UM- 



lated with THERMINATOR [231. With the freeze-out hyper- 
surfaces generated in this work, the collision rate after 
freeze-out is not very large. For the obtained hypersur- 
faces the number of pionic trajectory crossings at the 
distance corresponding to the cross section for the pion 
collisions is about 1.5-1.7. Hence the single-freeze-out 
scenario 



24| seems to be a fairly good approximation for 
the present case. The use of hadronic afterburners for 
elastic collisions has been described, e.g., in [25l l26l[27j. 

Our results for central and mid-peripheral collisions 
with the centrality classes adjusted to the available 
PHENIX [13, m and STAR ^ data are shown in Figs.[T] 
and [2] (darker lines/bands). We note a uniform agree- 
ment for all soft phenomena studied. In particular, the 
transverse-momentum spectra, the pionic elliptic-flow, 
and the HBT radii, including the ratio -Rout/^side, are 
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described within 10% or better. The HBT results for 
peripheral collisions are of similar quality as in Fig. [21 

The discussed results have been obtained with an early 
start of hydrodynamics, at the proper time r — 0.25 fm. 
It is unlikely that the system should equilibrate so early. 
Next, we argue that one may delay the starting point of 
hydrodynamics to realistic times by the inclusion of the 
partonic free-streaming between the initial proper time 
To = 0.25 fm when the partons are formed and some later 
time, T, when hydrodynamics starts. Similar ideas have 
been described by Sinyukov et al. in Refs. [H,!!^. Thus 
the global picture is as follows: early phase (CGC) gener- 
ating partons at time tq - partonic free streaming until r 
- hydrodynamic evolution until freeze-out at temperature 
Tf — free streaming of hadrons and decay of resonances. 

Massless partons are formed at the initial proper time 



^0 = 

speed of light until the proper time when free streaming 
ends, T — VP 
pidities rjQ - 
kinematics 



and move along straight lines at the 



^2 _ 2;2. We introduce the space-time ra- 
i log and 77 = i log . Elementary 



to+zo ' 2 t+z ■ 

[285 links the positions of a parton on the 
initial and final hypersurfaces and its four-momentum 
pf^ = {prCOshY, pt cos (j), px sin cj), pTs'mhY) , where Y and 
Pt are the parton's rapidity and transverse momentum: 



Tsinh(77 — Y) — Tosinh(77o — Y), 
X = xo + A cos 0, y = yo + A. sin 

A = ^ 
coshy 



(2) 



TCOsh(F — f]) — \l Tq + r2sinh^(F — ry). 



Thus the phase-space density of partons at the proper 
times To and t are related, 



dY (PpTd-qdxdy 



dr]Qdxodyo 



d^A^(To) 



dY d'^PTdrjodxodyo 



6{i]Q —Y — arcsinh[ — sinh(77 — Y)]) x 
To 

6(x — Xo — A cos (j))5{y ~ yo — A sin 0) 



(3) 



It is reasonable to assume a factorized boost-invariant 
form of the initial distribution of partons. 



d^N{To) 



dY d'^PTdrjQdxQdyci 



nixQ,yo)F{Y - t^cPt) 



(4) 



where n is the transverse density of Eq. ([T]). When 
the emission profile F is focused near Y — rjo, for in- 
stance ~ exp[-(y - r7o)^/(2(T^)], with ct 1, and if 
T ^ To, then the kinematic condition ^ transforms it 

into F ^ exp ^— arcsin 

so sharply peaked that effectively F 
Eq. (O yields 



^ s'm{Y — 1]) /(2a-^) j which is 
(5(y - 77). Then 



d^N{T) 



dYd'^pxdrjdxdy 



n(x — At cos (t>,y — At sin (jj) x 
5{Y~ri)fipT). (5) 



where At = t—tq and /(pt) is the transverse momentum 
distribution. The energy-momentum tensor at the proper 
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FIG. 3: (Color online) Sections of the energy-density profile e 
(gaussian-like curves) normalized to unity at the origin, and 
of the velocity profile v = \/Vx + Vy (curves starting at the 
origin), cut along the x axis (solid lines) and j/-axis (dashed 
lines) . The initial profile is from Eq. ([1} for centrality 20-40% 
at To — 025 fm. The e profiles are for t = tq = 0.25, 1, and 
2 fm, while the velocity profiles are for t = 1 and 2 fm, all 
from bottom to top. We note that the flow is azimuthally 
asymmetric and stronger along the x axis. 



time T, rapidity 77, and transverse position {x, y) is given 
by the formula 



T^"" = j dYd'pT 



d^N{i 



(6) 
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^0 

^ cosY? rj coshrj COB (f) coshr7sin0 cosh77sinh77 ^ 

cosh77 cos 4> cos^ 4> cos (j) sin cos 0sinh77 
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where A is a constant from the pt integration. Due to 
boost invariance the further calculation may be carried 
at 77 = 0. Next, we assume that at the proper time t 
the system equilibrates rapidly. We thus use the Landau 
matching condition, 



T^''{x,y)uy{x,y) = e{x,y)g^"'u^{x,y), 



(7) 



which states that at each point the four-velocity of the 
fluid, u^' = {l,Vj,,Vy,Q)/yJl - v'^, is such that T'^^ is di- 
agonal in the local rest frame of the matter. The position- 
dependent eigenvalue e is identified with the energy- 
density profile. The result of solving Eq. ([7]) with T''" 
from ^ for At = 0.75 fm and the initial profile ^ for 
c=:20-40% is shown in Fig. \3\ The curves show the sec- 
tions along the x and y axes of e and the velocity v at To 
(no free streaming) and t = 1 and 2 fm. Obviously, at 
To we find e{x, y) = 2Tr An{x, y). We note that, naturally, 
the profile spreads out as the time progresses. Impor- 
tantly, this effect is faster along the shorter This 
is clearly indicated by the velocity profiles along the x 
and y axes. Thus the flow generated by free streaming 
is azimuthally asymmetric. It can also be obtained upon 
the low At and low x, y expansion, where straightforward 
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4: (Color online) Study of azHBT: The quantities 
fiLn ^/flL^ nand i?out,2/^sidc,o from the model (bands) and 
experiment |30] (points), plotted as functions of the transverse 
momentum of the pion pair. 
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The results of the hydrodynamic calculation starting 
from free-streaming down to t = 1 fm and readjusted 
initial temperature at the proper time t, followed by 
THERMINATOR simulations, are shown in Figs. [1] and |1] 
with lighter curves/bands. We notice very similar results 
to the case with no free streaming, and again a proper 
description of the data. Larger free-streaming times spoil 
this agreement, as the flow becomes too strong. We re- 
mark, however, that the departure from the S{Y — rj) 



condition weakens the flow, which allows for larger val- 
ues of T. The basic conclusion here is that the partonic 
free streaming may be used to delay the start of hydrody- 
namics. The physical results are basically unaltered, as 
the dispersion of the density profile, resulting in milder 
hydrodynamic development of flow, is accompanied by 
the buildup of the initial flow. 

Having properly described the HBT radii, we may now 
deal "as a bonus" with the azimuthally sensitive HBT 
interferometry [s^] and consider the averages over the 
azimuthal angle, R^2{^t) — {Ri {kr , f) cos{2ip)) , where 
i — side or out (i?2,iong = 0). The results shown in Fig. S] 
display a remarkable agreement between our model and 
the data (the curves with and without free streaming 
overlap). 

In conclusion: 1) With a proper choice of the initial 
profile (e.g. Gaussian) one may obtain a uniform de- 
scription of soft observables (pr-spectra, V2, and pionic 
HBT radii) at RHIC with the boost-invariant inviscid hy- 
drodynamics with a realistic equation of state and thus 
solve one of the "RHIC puzzles" . 2) Partonic free stream- 
ing generates initial transverse and elliptic fiow. It may 
be used to delay the start of the hydrodynamic phase. 
3) Azimuthally-sensitive HBT is described in agreement 
with the data. 4) The complete treatment of resonances 
is important, as is well known from statistical models. 
5) When other effects (viscosity, departure from boost 
invariance) are incorporated, the important role of the 
initial condition should not be overlooked. 
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